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Hermann Feshbach predicted fifty years ago [l| that 
when two atomic nuclei are scattered within an open en- 
trance channel — the state observable at infinity, they 
may enter an intermediate closed channel — the locally 
bounded state of the nuclei. If the energy of a bound state 
of in the closed channel is fine-tuned to match the rela- 
tive kinetic energy, then the open channel and the closed 
channel "resonate" , so that the scattering length becomes 
divergent Q. We find that this so-called Feshbach reso- 
nance phenomenon not only exists during the collisions of 
massive particles, but also emerges during the coherent 
transport of massless particles, that is, photons confined 
in the coupled resonator arrays We implement the 
open and the closed channels inside a pair of such arrays, 
linked by a separated cavity or a tunable qubit. When 
a single photon is bounded inside the closed channel by 
setting the relevant physical parameters appropriately, 
the vanishing transmission appears to display this pho- 
tonic Feshbach resonance. The general construction can 
be implemented through various experimentally feasible 
solid state systems, such as the couple defected cavities 
in photonic crystals. The numerical simulation based 
on finite-different time-domain(FDTD) method confirms 
our conceive about physical implementation. 

The phenomenon of Feshbach resonance has been 
found in many physical systems over the years, such 
as the electron scattering of atoms 0] and diatomic 
molecules More recently, the development of laser 
cooling technologies has enabled the observation of low- 
energy Feshbach resonance in ultra-cold atoms @, 0, 0] 
and Bose-Einstein condensates (BEC) [^, [l(|. These 
experiments have helped verify the simulation of vari- 
ous theoretical predictions of condensing phenomena in 
solid state systems 111 ]. The latter in particular exem- 
plifies the resonance phenomenon as a means for adjust- 
ing inter-atomic coupling in reali zing various quantum 
phases ranging from BEC to BCS |l2l [l3j . 

On the other hand, since their first discovery by von 
Neumann and Wigner F| , bound states have been stud- 
ied in a general continuum [l5| and the emergence of a 
bounded energy level has been verified by various mod- 



els [16|, [13, 113], extracting from the simple coupling of 
a discrete level with a continuum. Quasi-bound states 
have been predicted in tight-binding fermionic quantum 
wires [ljj for localized fermions and in the optical cou- 
pled resonator arrays 0, 0, 21 1 for confined photons. 
It has also found various applications in many quantum 
optical devices 22 , [23|, [24| , including the single-photon 



transistor [25|. Under this retrospect, we ask if it is possi- 
ble to control exactly when the photons become bounded 
and unbounded through an external parameter, viz. to 
implement an optical version of Feshbach resonance. 

The desired resonance between the bound and the un- 
bound states can be found inside a pair of parallelly 
placed coupled resonator arrays 0], a series of consec- 
utively placed optical microcavities that entrap photons 
and allow photon-hopping from one of the cavities to its 
closest neighbors at left and right. The two arrays are 
connected by a central cavity that acts as a quantum 
controller and couples separately to one cavity in each 
of the arrays, as shown in Fig. [TJa), forming an H-shape 
system. We designate the upper array as array A with 
the Hamiltonian 

H = u A ^ a ] a j + (Ja ^2 a ] a J+i + 9Aa c^ + h.c), (1) 

3 j 

where the second summation describes the tight-binding 
hopping of photons between neighboring cavities with 
hopping coefficient Ja while the first one accounts for 
the static photon occupations in the cavities. a 3 - denotes 
the annihilation operator of the bosonic mode for the j- 
th cavity field and we assume that the mode frequency 
lua of each cavity field is identical. The last term is the 
interaction between the central cavity C of single mode 
frequency O and the zeroth resonator of array A with cou- 
pling strength gA- c is corresponding annihilation bosonic 
operator for central cavity. The lower array is designated 
as array B, the Hamiltonian Hb for which is no differ- 
ent from Eq. fl} except the change of the bosonic mode 
operator to bj ,the hopping coefficient to Jb , the mode fre- 
quency to lub and the coupling strength to <?b- Then, the 
total model Hamiltonian reads H = Ha + Hb + fic^c. 

One experimental available system of the model is 
based on photonic crystal, which will be described in de- 
tails below. We have to point out that there is a draw- 
back in our setup: the inter-cavity coupling cannot be 
externally manipulated, once the photonic-crystal based 
metamaterial is fabricated as an all-optical chip. How- 
ever, for single photon transferring, the role of central 
cavity as that of a qubit with two levels |e) = |1) and 
\g) = |0), corresponding to the single photon state and 
vacuum of the central cavity. This identification moti- 
vates us to use a qubit controller replacing the central 
cavity equivalently in the single photon case. Therefore, 
the phenomenon predicted here can be realized in the hy- 
brid system of photonic crystal and quantum dots with 
external controllable parameter. Our general construe- 
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FIG. 1: The photonic Feshbach resonance based on the cou- 
pled resonator arrays: (a) schematic of the H-shape system 
consisting of two coupled resonator arrays, array A and array 
B, connected by a qubit controller of level spacing Q. Each 
of the cavity in the two arrays is characterized by the cavity 
field mode it confines and hereby indicated by the bosonic 
operator aj or bj and indexed by the relative distance j from 
the zeroth cavity with which the qubit controller couples; (b) 
probability distribution of a single photon in a coupled res- 
onator array for a quasi-bound state. When the level spacing 
SI of the qubit controller matches certain values, the probabil- 
ity amplitude would vanish at the two ends of the array; (c) 
the energy state distribution for array a, a £ {A, B}, relative 
to the common eigenfrequency iu a of the cavity fields: a pair 
of a continuum and a discrete level above and below iv a . 



tion can also be implemented through the circuit QED 
system [2o| including two coupled superconducting trans- 
mission line cavity linked by charge or flux qbit. The 
replacement of the central cavity by a controllable qubit 
can overcome the drawback mentioned above. 

In the following , we use the Jaynes-Cummings cou- 
plings g a c^ \g) (e| + h.c (a = A, B) to modeling the 
central cavity couplings in the single photon case. 
The HUbert space is spanned by the tensor product 
(l e ) j Iff)} ®j l n Aj) ®j \ n B,j) where \n a j) denotes the 
state of array a with its j-th cavity being occupied with 
n a ,j number of photons. When separated from the other 
and studied individually, each coupled resonator array 
is described by the subsystem Hamiltonian H a and pos- 
sesses two bound states, reminiscent that of the Feshbach 
resonance. The states are the particular superposition of 
eigenstates for the Hamiltonian H a , comprised by a sub- 
set of the basis vectors described above. For array A, the 
state is namely \ip A ) = J2 3 u A,g{j) Iff) lj, 0) +u A , e \e, 0, 0) 
where only the excited state of the central cavity and a 
single-photon excitation in one of the cavities, as indi- 
cated by the lj symbol, are included. The coefficients 
UA,g(j) and UA, e in the equation constitute the spectrum 
of probability distributions of these states. That of array 



B takes a similar form with amplitudes UB, g (j) and «B,e- 
To see whether there are bounded single-photon states 
within their individual coupled resonator array, we can 
solve the discrete-coordinate scattering equation associ- 
ated with their corresponding eigenvalue E for the prob- 
ability spectrum, 

[E -(j a -V a (E)}u a ,g(j) = -J a [u at g(j + l)+u at g(j -1)] 

(2) 

where the term V a (E) = g^8jo/(E — O) on the left hand 
side is contributed by the JC type interaction between the 
central cavity and the coupled resonator array. V a (E) is 
a resonate potential that depends on the eigenenergy E. 
In the continuous limit of the coordinate j — > x, this term 
reduces to a S -type potential 
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The 5-type potential forms a confining barrier to the 
transportation of single photon in the coupled resonator 
array and informs a bounded single photon within, simi- 
lar to those in the models proposed in Refs. 2(3, 21 1. 

It has a singularity at E being equal to the level spac- 
ing fi, leading to a quasi-plane- wave type solution Q 



to Eq. 



/(j) = Ca exp(-zK Q |j|) where C a de- 



notes a constant and the wave number is complex, n a = 
K>a,R~ ii^a,i- The imaginary part K a j of the wave number 
can admit a positive value and for the non-zero coupling 
g a , resulting in a decay of the probability distribution of 
single-photon states over the discrete spatial coordinate 
j. The vanishing probability amplitude towards the ends 
of the arrays, i.e. along with \j\ oo, demonstrates the 
existence of a bound state of a single photon, as shown in 
Fig. [ljb). For this system, continuum band has a band- 
width of 4J a and their paired discrete levels, denoted 
respectively by E a + and E a -, are gapped from cither 
below or above, as illustrated in Fig. [TJc). Reverted to 
the conventional language of atomic scattering, the con- 
tinua of cigcncncrgies can be considered open channels of 
multiple admissible energy states in the continuous range 
uj a — 2J a < E < uj a + 2J a . Out of this range, the energy 
states can only admit two discrete levels that associate 
with a non-real k a , representing closed channels or bound 
states. These two discrete levels 



E a ± = r2 i 



Oo 



4J2 



(4) 



are exactly solved from the above discrete-coordinate 
scattering equation. The dependence of these two bound 
state energies on the various system parameters gives 
hints to their potential of tunability and controllability. 
The readers familiar with the approaches by Lee [lq . 
Fano [It]], and Anderson [l8[ shall also find our result 
here familiar. 

The next logical step is to study how the reso- 
nance phenomenon arises when two individual arrays 
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are paired to form the H-shape system. The scattering 
state \<p) of the H-system for single-photon reads \ip) = 



E,' [ u A,gU) Iff. ij. ) + u B,gU) Iff. , h)} 



,0,0). The 



probability amplitudes UA,g{j) and Ub,s(j) °f the pho- 
tonic occupation among the cavities and u e of the atomic 
excitation in the system can still be analyzed through 
the time- independent Schrodinger equations, leading to 
a pair of algebraic scattering equations similar to Eq. (|2j). 
one for the probability amplitudes in each array. The 
distinction of the case here lies in the adding W a (E) = 
gagaU a , g {0))/{E — SI) in the right hand side with a in- 
dexing either A or B. We have used a to indicate the dual 
array relative to a; namely, when a indexes A, then a 
indexes B and vice versa. This term reflects a potential 
again contributed by the interaction of each array with 
the central cavity. However, because of the coupling be- 
tween the central cavity and the dual array, additional 
contribution from the dual array has to be considered. 

The set of solutions to the pair of scattering equations 
are many. The portion we are concerned with are those 
illustrating a simultaneously existent set of open channel 
and closed channel in the two coupled resonator arrays. 
We select one of the particular cases when a single photon 
is inserted into an open channel in array A from the left 
for this purpose. The distribution of this single photon 
in the array is then described by a plane wave, ua(j) = 
exp(z/cj)-|-rexp(— ikj) (j < 0); sexp(ifcj) (j > 0). s and 
r denote, respectively, the transmission and the reflec- 
tion coefficients of the optical plane wave, indicating the 
scattering of photon by the effective potential Eq. {3J at 
the zeroth resonator in the one-dimensional coordinate 
space. Meanwhile, the distribution amplitude for the 
single photon in array B can be quasi-plane-wave type, 
u b(j) = Cb exp {— iK\j\}, with a complex wave number 
k, and indicate a closed channel, same as that of the 
individually discussed case. These two distributions in 
the paired arrays, when combined through the coupled 
scattering equations, give rise to a unified dual-channel 
coupling equation 

(1 - s) sin fc— = C B sin re— = 9 ^ S ^ 9B ^ . (5) 

Therefore, the optical dual-channel resonance occurs 
when there exists a solution of real k and complex re to 
Eq. ([5]) and the eigenenergy E of the photon in array A 
matches either of the discrete energy levels Eb± of array 
B. The process is illustrated in Fig. [2] for two particu- 
lar cases with E matching Eb+ in Fig. [2] (a) or Eb- in 
Fig. EJb). The two possibilities of channel resonances is 
further illustrated in Fig. [2] (c) with the equivalent poten- 
tial of array A as a function of the resonator position j. 
The zeroth resonator locates the position of local minimal 
energy for both the open channels and the closed chan- 
nel, reflecting the potential barrier set up by the central 
cavity. The dual-channel resonance occurs as well when 
the roles of array A and array B arc exchanged. 
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FIG. 2: Diagrams of the energy state distributions of the 
single photon in array A and array B, showing the process 
of Feshbach resonance between an open channel and a closed 
channel. Two particular resonance cases exist for an incident 
photon inserted into array A: (a) the photon energy level in 
the continuum band in array A is resonant with the upper 
discrete level Eb+ in array B; (b) the photon energy level in 
the continuum band in array A is resonant with the lower 
discrete level Eb- in array B. (c) A profile plot illustrating 
the same two cases of channel resonances relative to the cavity 
position j the photon occupies. 



The criteria of the dual-channel resonance can be met 
if the transmission coefficient s in Eq. © vanishes. This 
condition implies the circumstance where the incident 
photon in array A is totally reflected or scattered by 
the potential barrier set up by the central cavity at 
position j = 0. In other words, the level spacing SI 
of the controller becomes our tuning parameter for the 
photonic Feshbach resonance. Written in terms of the 
other variables in the dual-channel coupling equation s = 
—2iCbJb{E — fi) sin re/ (<7a3b) — Cb.9b/.9a, the transmis- 
sion coefficient vanishes when E = SI — gj|/(2i Jb sin re), 
leading to a complex wave number re as expected. The 
complete reflection in the open channel can be under- 
stood as the divergence effect of s-wave scattering length 
for the usual Feshbach resonances in three-dimensional 
space reduced to a version in one-dimensional space. 
Eliminating the various variables, the transmission co- 
efficient can be expressed as a function of the incident 
energy E, i.e., 
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s = 



\f a (E)[F a {E)-G-(E)] \ £>w B t^b 



[F A (E)[F A (E)-G+(E)} \ E<uj b 



2J B 
2J B 



where we have 



^(E-u a y- 

glF B (E)/((E 
flection coefficient 



used 
Ut with 



the 



shorthands F a (E) = 
{A, B} and G ± {E) = 
Q)F B (E) ± g' B ). The norm-squared re- 



a e 

,2 s 



|1 — s\ is plotted against the photon 



energy in families of varying level spacing £1 and coupling 
constant g A of the qubit controller in Fig. [3] The pho- 
ton encounters two kinds of characteristic points while 
propagating through array A. The first one is an indiffcr- 
cntiablc turning point where s = 1 or E = lu b + 2 J B . The 
potential barrier becomes transparent and the photon is 
completely transmitted because of the matching coupling 
between the qubit controller and the dual array B. The 
second one is the maximum point where the photon is 
fully reflected when the transmission coefficient is van- 
ishing s = 0. We hence see the shifting of this peak while 
f2 is varied. The reliance on the coupling coefficient g A 
determines the width of the peaking. 

Next we numerically examine the feasibility of our the- 
oretical prediction on a two-dimensional photonic crys- 
tal 27|, [28j. The crystal is made up of a square lat- 
tice of high-index dielectric rods of radius 0.2a, 0.1a 
and 0.05a, where a is the lattice spacing. The artifi- 
cial design is made by two parallel waveguides of cou- 
pled defected cavity arrays linked through a central de- 
fected cavity on the two-dimensional photonic crystal, as 
illustrated in Fig. [4ja). The two resonator arrays [29[ 
is constructed with different frequencies, inter- resonator 
tunneling rates, and coupling strengths with the central 
cavity. 

For this photonic crystal, the material of all the rods 
is assumed to be silicon, with a dielectric constant e = 
11.56, and the background is filled by air. We make 
the simulation of the designed structure with the finite- 
difference time-domain (FDTD) method [3(| in freely 
available Meep code [3l[ . The steady field vector of the 
incident wave at frequency luq = 0.3628 x 27rc/a is plot- 
ted in Fig. BJb), with red showing positive amplitudes 
pointing out from the plane and blue negative amplitudes 
into the plane. The wave travels horizontally from left to 
right and hence, according to the convention for charac- 
terizing photonic crystals, carries a transverse-magnetic 
(TM) polarization. The notice-worthy region is located 
at the center where the highly-saturated colors indicate 
a localized bounded photon from the lower waveguide. 
Moreover, the blank portion in the upper waveguide in- 
dicates a completely reflected wave. Finally, we point out 
that, though the numerical simulation based on FDTD is 
of classical, but the weak light calculation can also reflect 
the single photon nature with the intensity distribution 
illustrated in Fig. Hlb) , which is only relevant to the first 
order coherence function. 
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FIG. 3: Plots of the norm-squared reflection coefficient |1 — s| 2 
against the eigenenergy E of a propagating photon in the 
coupled resonator array A of the process illustrated in Fig. 
[2^a). Two tuning parameters are varied: (a) the level spacing 
SI of the qubit controller; and (b) the coupling coefficient g\ 
between the qubit controller and the zeroth cavity in array 
A. Other parameters are chosen to be: Ja = 1, Jb = 0.5, 
uja = 2, ujb = 1 and <7b = 0.7. The incident energy ranges 
from to 4 and the continuum band for the array B is set 
[0, 2]. The corresponding bound state energies are marked as 
Ebi,Eb2,Eb3 in (a) and E b in (b). 



In conclusion, we have shown the existence of a pho- 
tonic bound state in a qubit-controllcd coupled res- 
onator array and predicted photonic Fcshbach resonance 
emerges from a pair of these coupled resonator arrays 
coupled in an H-shape fashion. An FDTD simulation of 
the system implemented on a photonic crystal has ver- 
ified the validity of the proposal. The resonance phe- 
nomenon arises from the dual-channel coupling between 
an unbound state in one array and a bound state in the 
other, the occurring moment of which is indicated by a 
total reflection of an incident photon in the array. Our 
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FIG. 4: An experimental protocol based on a photonic crystal 
made up of silicon rods of radius 0.2a. (a) The structure of 
the design: the upper waveguide is implemented by removing 
a row of original rods and substituting with a set of rods of 
radius 0.1a and spacing dup = 4a. The lower waveguide, D = 
la apart from the upper one, is constructed by removing three 
rods out of every five rods, i.e. lattice spacing doowN = 5a. 
The central cavity is created by reducing the radius of three 
vertically-placed rods between the two waveguides to 0.05a. 
(b) Plotting the steady electric field vector for an incident 
wave of frequency luq — 0.3628 x 2nc/a with TM-polarization. 



analysis for the resonant scattering process was carried 
out for the single photon case and did not rely on the 
photonic statistics. Our prediction here is thus applica- 
ble to fermionic models, such as the electron transporta- 
tion along a H-shape array of quantum dots. For the 
case where multiple photons are assumed to exist in the 
arrays, Bethe-ansatz must be used for the analysis and 
we shall defer its discussion in a future work. 

C.P.S. acknowledges the helpful discussion with S. 
Yang, Peng Zhang and X. H. Wang. This work 
is supported by NSFC No. 10474104, No. 60433050, 
and No. 10704023, NFRPC No. 2006CB921205 and 
2005CB724508. 



[1] 
[2] 

[3] 



Electronic address: dhui@itp.ac.cn 

Electronic address: suncp@i tp.ac.cn| 

Feshbach, H. Unified theory of nuclear reactions. Ann. 

Phys. (N.Y.) 5, 357-390 (1958). 

Pethick, C. J. & Smith, H. Bose-Einstein Condensation 
in Dilute Gases (Cambridge University Press, 2002). 
Zhou, L., Gong, Z. R., Liu, Y.-X., Sun, C. P. & Nori, F. 
Controllable scattering of a single photon inside a one- 
dimensional resonator waveguide. Phys. Rev. Lett. 101, 
100501 (2008). 



K 
[5, 

ir 

[8, 
[9 

[10 

[11 
[12 

[13 

[14 

[is; 

[16 

[17; 

[18 
[19 

[20 

[21 

[22 
[23 
[24 



Schulz, G. J. Resonances in electron impact on atoms. 
Rev. Mod. Phys. 45, 378-422 (1973). 
ibid. Resonances in electron impact on diatomic 
molecules. Rev. Mod. Phys. 45, 423-486 (1973). 
Courteille, P., Freeland, R. S., Heinzen, D. J., van Abee- 
len, F. A. & Verhaar, B. J. Observation of a feshbach 
resonance in cold atom scattering. Phys. Rev. Lett. 81, 
69-72 (1998). 

Roberts, J. L., Claussen, N. R., Burke, J. P. Jr., Greene, 

C, Cornell, E. A. & Wieman C. E. Resonant magnetic 
field control of elastic scattering in cold 85Rb. Phys. Rev. 
Lett. 81, 5109-5112 (1998). 

Vuletic, V., Kerman, A. J., Chin, C. & Chu, S. Obser- 
vation of low-field feshbach resonances in collisions of ce- 
sium atoms. Phys. Rev. Lett. 82, 1406-1409 (1999). 
Inouye, S., Andrews, M. R., Stenger, J., Miesner, H. J., 
Stamper-Kurn, D. M. & Ketterle, W. Observation of Fes- 
hbach resonances in a Bose-Einstein condensate. Nature 
392, 151-154 (1998). 

Timmermans, E., Tommasini, P., Hussein, M. & Ker- 
man, A. Feshbach resonances in atomic Bose-Einstein 
condensates. Phys. Rep. 315, 199-230 (1999). 
Jin, D. S. &i Regal, C. A. Proceedings of the Interna- 
tional school of physics Enrico Fermi, course CLXIV (IOS 
Press, Amsterdam, 2008). 

Loftus, T., Regal, C. A., Ticknor, C, Bohn, J. L. & Jin, 

D. S. Resonant control of elastic collisions in an optically 
trapped fermi gas of atoms. Phys. Rev. Lett. 88, 173201 
(2002). 

Zirbel, J. J., Ni, K.-K., Ospelkaus, S., D'Incao, J. P., 
Wieman, C. E., Ye, J. & Jin, D. S. Collisional stability 
of fermionic feshbach molecules. Phys. Rev. Lett. 100, 
143201 (2008). 

von Neumann, J. & Wigner, E. Uer merkwiirdige diskrete 

Eigenwerte. Z. Phys. 30, 465-467 (1929). 

Friedrich, H. & Wintgen, D. Physical realization of bound 

states in the continuum. Phys. Rev. A 31, 3964-3966 

(1985). 

Lee, T. D. Some special examples in renormalizable field 
theory. Phys. Rev. 95, 1329-1334 (1954). 
Fano, U. Effects of configuration interaction on intensi- 
ties and phase shifts. Phys. Rev. 124, 1866-1878 (1961). 
Anderson, P. W. Localized magnetic states in metals. 
Phys. Rev. 124, 41-53 (1961). 

Nakamura, H., Hatano, N., Garmon, S. & Petrosky, 
T. Quasibound states in the continuum in a two chan- 
nel quantum wire with an adatom. Phys. Rev. Lett. 99, 
210404 (2007). 

Gong, Z. R., Ian, H., Zhou, L. & Sun, C. P. Controlling 
quasibound states in ID continuum through electromag- 
netic induced transparency mechanism. arXiv:0805.3042 
(2008). 

Dong, H., Gong, Z. R., Ian, H., Zhou, L. & Sun, C. P. 
Intrinsic cavity QED and emergent quasi-normal modes 
for single photon. larXiv:0805 .3085 (2008). 
Shen, J. T. & Fan, S. Coherent single photon transport in 
a one-dimensional waveguide coupled with superconduct- 
ing quantum bits. Phys. Rev. Lett. 95, 213001 (2005). 
ibid. Strongly correlated two-photon transport in a one- 
dimensional waveguide coupled to a two- level system. 
Phys. Rev. Lett. 98, 153003 (2007). 

ibid. Coherent photon transport from spontaneous emis- 
sion in one-dimensional waveguides. Opt. Lett. 30, 2001- 



6 



2003 (2005). 

[25] Chang, D. E., S0rensen, A. S., Demler, E. A. & Lukin, 
M. D. A single-photon transistor using nanoscale surface 
plasmons. Nature Phys. 3, 807-812 (2007). 

[26] Wallraff, A., Schuster, D. I., Blais, A., Frunzio, L., 
Huang, R.-S., Majer, J., Kumar, S., Girvin, S. M. & 
Schoelkopf, R. J. Strong coupling of a single photon to a 
superconducting qubit using circuit quantum electrody- 
namics. Nature 431, 162-167 (2004). 

[27] Fan, S., Sharp asymmetric line shapes in side-coupled 
waveguide-cavity systems. App. Phys. Lett. 80, 908-910 
(2002). 

[28] Joannopoulus, J. D., Johnson S. G., Winn, J. N., & 
Meade, R. D., Photonic Crystal: Molding the flow of light 



(Princeton University Press, 2007). 

[29] Yariv, A., Xu, Y., Lee, R. K. & Scherer, A., Coupled- 
resonator optical waveguide: a proposal and analysis. 
Opt. Lett. 24, 711-713 (1999). 

[30] Tafiove, A., & Hagness, S. C, Computational Electro- 
dynamics: The Finite-Difference Time-Domain Method 
(Artech: Norwood, MA, 2000). 

[31] Farjadpour, A., Roundy, D., Rodriguez, A., Ibanescu, 
M., Bermel, P., Joannopoulos, J. D., Johnson, S. G., 
& Burr, G. Improving accuracy by subpixel smooth- 
ing in FDTD. Optics Letters 31, 2972-2974 (2006). 
|http://ab-initio. mit.edu/wiki/index.ph p/Meep 



